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Abstract. An improved forced- precession model of 
Comet 46P/Wirtanen is presented. The nongravitational 
motion of the comet has been investigated based on the 
forced precession model of the rotating cometary nucleus. 
The least squares method applied to observational equa- 
tions allows us to determine six basic nongravitational pa- 
rameters: A, 77, I, 4>, f p , and s together with six orbital 
elements. The solutions were obtained with additional as- 
sumptions: 1. The cometary activity with respect to the 
perihelion is asymmetric and, 2a. The time shift parame- 
ter which describes the displacement of maximum activity 
with respect to the perihelion could change its value dur- 
ing the investigated apparitions or 2b. The real activity of 
the comet changed during the time interval considered . 
All eight recorded apparitions of the comet were linked us- 
ing all astrometric observations covering the period 1948- 
1997 with a mean RMS residual of l'/6. According to the 
best solution, the nucleus of Comet Wirtanen is oblate 
along the spin-axis with a ratio of equatorial to polar ra- 
dius of about R a /Rb = 1.1. This precession model yields 
a value of 4.9 hrs/km for the ratio of P m t/R a - Assuming 
that the nucleus radius of 46P /Wirtanen most likely is in 
the range of 0.5-2.0 km we obtain a range 2.5 - 10.0 hours 
for the rotational period. 

Key words: Comets: individual: 46P/Wirtanen 



1. Introduction 

The short-period comet 46P/Wirtanen will be investi- 
gated for two years from 2011 during the ROSETTA 
cometary rendez-vous mission. Therefore, precise predic- 
tions of the future perihelion passages and space positions 
of the comet are crucial for successful planning of the mis- 
sion. 

Comet 46P/Wirtanen has been observed at eight re- 
turns during the almost hfty-year interval since its discov- 
ery in 1948 by C.A. Wirtanen. Two close approaches of the 



comet to Jupiter in 1972 and 1984 are responsible for sig- 
nificant changes of orbital elements. Apart from the plan- 
etary perturbations the comet is subject to nongravita- 



tional forces (Marsdcn & Williams 1997). The method of 
determining the nongravitational effects in the orbital mo- 



tion of the comet was proposed by Marsden et al.| (1973). 



In Marsden's formalism the three orbital components of a 
nongravitational force acting on the comet have the form: 



Fi = Aig(r), Ai — const for i = 1,2, 3, 



(1) 



where F±, F2, F$ represent the radial, transverse and nor- 
mal component of the nongravitational force, respectively. 
The analytical function g(r) simulates the ice sublimation 
rate as a function of the heliocentric distance r. 

In a previous paper ( jKrolikowska fc Sitarski 1996, 
hereafter Paper I) we concluded that the forced preces- 
sion model of the rotating nonspherical cometary nucleus 
adequately explains the variations of the nongravitational 
effects observed in the periodic comet Wirtanen and is 
suitable for making predictions of the future returns. In 
the present paper, we continue the study of nongravita- 
tional motion of the 46P/ Wirtanen including the astro- 
metric observations from the last apparition. 

A forced precession of the spin axis is caused by the 
torque which arises when the vector of the jet force (ex- 
erted by the outgassing on the nucleus) does not pass 
through the mass center of the nonspherical nucleus. The 
precession rate is a function of the nucleus orientation, 
the lag angle r\ of the maximum outgassing behind the 
subsolar meridian, the modulus of the reactive force, the 
nucleus oblateness s and the precession factor f p which de- 



pends on rotation period and nucleus size (Sekanina 1984 



1988). The orientation of the nucleus is defined by two an- 
gles: the obliquity / of the orbital plane to the cometary 
equator and the longitude <fi of the Sun at perihelion. The 
nongravitational force components (which vary along the 
orbit) are related to the angular parameters 77, /, 4> of the 
rotating nucleus by: 

F i (t)=A-C i (ri,I(t),v(t)+<p(t))-g(r), for i= 1,2,3, 
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where A = + F| + F%) 1 / 2 / g(r), v is the true anomaly 
of the comet and Ci denote the direction cosines which 
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Table 1. Orbital elements and physical parameters of the nucleus for forced precession models linking all apparitions 
of 46P/Wirtanen. Angular elements u>, f2, i refer to Equinox J2000.0. Parameters A, A 1 and A 11 are in units of 
10 -8 AU/day 2 , the precession factor f p is in units of 10 6 day/AU, and time shifts r, t\ and T2 are in days. 



Model 1 



Model 2a 



Model 2b 



Equinox: 
Epoch: 

T 

9 
e 
us 
Q 



A 

A 1 

A 11 

V 
/<> 

U 

s 

T 
Tl 

T2 

Res 



1997 03 14.15654 
1.06342371 
0.65681460 
356? 32749 
82? 20129 
11? 72345 



Orbital elements 



J2000.0 
1998 01 27 
ET = JD 2450840.5 
1997 03 14.15568 
1.06342736 
0.65680686 
356? 33091 
82? 19688 
11? 72330 



1997 03 14.15611 
1.06343483 
0.65680623 
356? 33058 
82? 19803 
11? 72321 



Nongravitational and 
nuclear parameters 



+0.62615 ±0.00523 



21? 02 ± 0?42 
136? 00 ± 1?13 
259? 60 ± 3? 82 
0.62528 ± 0.02536 
0.33419 ±0.01318 



1'.'82 



-0.76153 ± 0.00986 



15?90 ± 0?58 
154? 64 ± 0?52 
23? 00 ± 8? 98 
+1.7461 ±0.1242 
+0.12148 ± 0.04720 



-31.623 ± 1.691 
-11.236 ±0.529 



+0.80188 ± 0.00873 
+0.67797 ± 0.00722 
15? 48 ±0?66 
140? 83 ± 1?28 
318?32 ± 1?82 
±1.2337 ± 0.0694 
±0.10190 ±0.03421 

-23.476 ± 1.247 



1'.'59 



are time-dependent due the orbital and precessional mo- 
tion of the comet. We adopt Marsden's water production 
curve, g(r), because of its analytical form which is useful 
for orbital computations. Expressions for time variations 
of the angles / and (f>, the precession rate, and the di- 
rection cosines of the nongravitational force are given by 
Krolikowska et al.| (1998). 



The values of five precessional parameters: A, ij, I , 
00: fp were derived in Paper I; the sixth - oblateness of 
the comet nucleus - could not be determined because of 
poor observational material which then consisted of 67 po- 
sitional observations covering the period 1948-1991. Since 
that analysis, the new, high-quality data have become 
available from the latest apparition 1995-1997. This al- 
lows us to refine the previous forced precession model. 
The aim of this paper is to obtain the best model of the 
comet by linking all the astrometric observations made 
during the eight recorded appearances of the comet. 



2. Observational Material 

The distribution of the 67 observations of Comet 
46P/Wirtanen carried out before 1995 is given in Table 1 
of the Paper I. At present, 247 positional observations (494 
residuals) are available from the recent apparition cover- 
ing the time interval of 26 June 1995 - 30 December 1997. 
Since the observations are distributed highly nonuniformly 
over the apparitions, we replaced more than two observa- 
tions of the same day by one average comet position (so- 
called normal place) , which yields 300 residuals for the last 
return of the comet. Combining mean residuals 1'.'96, 1'.'67, 
l'/52, 2'.'29, 2'.'44, T/97 and l'.'Ol calculated separately for 
apparitions in 1948, 1954, 1961, 1967-1975 (two appari- 
tions), 1986, 1991 and 1996, respectively, w e obtained an 
a priori mean residual ( Bielicki fc Sitarski 1991) of l'/38. 
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3. Forced precession model of the rotating 
nonspherical cometary nucleus 

The equations of the comet motion have been in tegrated 



numerically by the recurrent power series method (Sitarski 



1979, 1984), taking into account all the planetary pertur- 
bations. A detailed description of the forced precession 
model is given in Paper I. The first set of equations con- 
tained six unknown parameters A, ij, In, 4>q, f P and s 
which were determined simultaneously with six orbital el- 
ements from the observational equations by an iterative 
least squares process. This solution is given in the first 
column of Table 1 (hereafter Model I). The positive value 
of s implies that - according to our model - 46P / Wirtanen 
rotates on its shorter axis. 

However, a comparison of the mean residual of this solu- 
tion with the a priori mean residual shows that a better 
solution is required for the motion of the comet. There- 
fore, we no longer make the assumption that the activity 
reached a maximum exactly at perihelion. We have found 
(from our orbital calculations) that between the appari- 
tions of 1986 and of 1991 (one revolution after the sec- 
ond closest approach of the comet to Jupiter) something 
changed in the cometary activity and/or in the geome- 
try of the sublimation area on the nucleus surface. These 
changes can be modelled by setting the discontinuity of 
parameters A and/or r close to the moment of aphelion 
passage between 1986 and 1991: 



.4 



A 1 for t < 1989.0 
A 11 for t > 1989.0 



and/or r 



where parameter r represents a time shift of the maximum 
value of g{r) with respect to the perihelion time. Thus, the 
function g(r) is then replaced by g(r'), where r' — r(t — r) 
( f5ekanina| 1988, |Sitarski| 1994). 

However, it turned out that due to the relatively poor 
observational material (especially that collected before 
1995), only two additional parameters could be numeri- 
cally determined. Thus, we were limited to introduce the 
discontinuity of r and A in two separate models, denoted 
as Model 2a and Model 2b respectively in Table 1. 
From this reason, the final two solutions are described by 
eight parameters: 

A, rj, I , <t>o, f P , s, n t 2 (Model 2a), 
and 

A 1 , A", 77, I , fa, f P , a, t (Model 2b), 
which were determined from the observational equations 
along with the orbital elements. These solutions, charac- 
terized by mean residuals of l'/68 and l'/59, are given in 
Table 1. Both models give a prediction of perihelion pas- 
sage in February 2008 (the last before the ROSETTA 
rendez-vous withthe comet) with a time dispersion of 
0.038 day. 

The motion of the nucleus rotation axis represented 
by angles / and 4> an( i the qualitative variations of the 
nongravitational force components Fi, F3, acting on 




ri for t < 1989.0 
r 2 for t > 1989.0 



I960 1980 3000 

DATE (years) 



Fig. 1. Temporal variation of the angle I - a, <f> - b, and 
components Fi,F 2 , F3 of the nongravitational force F - 
c, for Comet 46P/ Wirtanen due to the spin-axis forced 
precession of the comet's nucleus. Solid curves represent 
Model 2b and dashed curves - Model 2a. Dotted horizontal 
lines in b indicate moments when the sunlit pole of nucleus 
in perihelion changes from 'northern' to 'southern' and 
vice versa. 



the comet during its successive returns to the Sun are 
presented in Fig. 1 for Model 2b (solid curves). Dashed 
curves in the same figure represent the variation of I and 
<f> for Model 2a; the variations of F 1 ,F 2 ,F 3 are undistin- 
guishable from those for Model 2b. This means that the 
same character of the nongravitational accelerations can 
be interpreted by two different models of cometary activ- 
ity. The time variations of /, <fi and F\,F 2 , F3 are a conse- 
quence of the forced precession and orbital changes caused 
by two close approaches to Jupiter in 1972 and 1984. These 
close encounters with Jupiter reduced the perihelion dis- 
tance from 1.63 AU in 1947 to 1.08 AU in 1986, which 
caused an increase of the nongravitational force during 
this period (Fig. lc). Both solutions imply a retrograde ro- 
tation of the cometary nucleus (I >90°). This is in agree- 
ment with the Jorda & Rickman (1995) conclusion from 
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Table 2. Evolution of orbits of 46P/Wirtanen according to Model 2b given in Table 1; perihelion distance q is in AU, 
period P is in years, angular parameters uj, f2, i referto Equinox J2000.0. 



T 



Epoch 



1997 03 14.14929 
2002 08 26.63663 
2008 02 02.13771 
2013 07 08.64672 



1.06376188 
1.05876883 
1.05750021 
1.05209806 



0.65673107 
0.65780736 
0.65807737 
0.65923863 



5.46 
5.44 
5.44 
5.42 



356? 34091 
356? 39605 
356? 34611 
356? 33829 



82? 20593 
82? 17493 
82? 17490 
82? 16315 



11? 72249 
11? 73837 
11? 73953 
11? 75731 



140? 18 
145? 09 
149? 17 
161?23 



309? 42 
357?62 
38? 25 
85? 46 



1997 03 13 
2002 09 03 
2008 01 15 
2013 07 07 



the positive lightcurve asymmetry combined with the neg- 
ative rate of comet orbital period change. 

4. Discussion 

The present precession models impose some constraints 
on the physical properties of the nucleus of Comet 
46P/Wirtanen. Using these solutions it is possible to cal- 
culate values of Sekanina's (1984) torque factor f tor — 
fp/s. Using this parameter we are able to calculate the 
value of the P ro t/Ra ratio from the formula: P ro t/Ra = 
Am ■ f P l (5s), where the equatorial radius, R a , is related to 
the oblateness and the polar radius, -R&: Rb/R a — 1 — s. 
Our models give: 



0.75 ± 0.01 for Model 1 
P rot /R a = \ 5.80 ± 1.85 for Model 2a 
4.88 ± 1.36 for Model 2b 



(2) 



where P rot is in hours and R a in kilometers. Tak- 
ing into account observations of the nucleus radius of 
46P/Wirtanen it is possible to verify our models. 
From several estimates of the nucleus radius of Comet 
Wirtancn it appears that this comet possesses one of the 
smallest nuclei detected so far. From photographic mag- 
nitudes Jorda & Rickman (1995) obtained a radius of the 
nucleus "of less than 2 km, probably 1.5-1.8 km". This 



is in agreement with Almeida et al. (1997), who derived 
a lower limit of the nucleus radius of 1.0 km from visual 
magnitudes of the comet. However, systematically smaller 
sizes of the nucleus are reported from CCD photometry. 
CCD observations done by Bochnhardt et al. (1997) give 



an upper limit of 0.8 km and Lamy et al. (1998) obtained 
a mean effective radius of 0.60 ± 0.02 km. Such small 
radii imply that a large fraction of the surface of the nu- 
cleus is active. The active area of 1.80 km 2 obtained by 



Farnham & Schleicher (1998) covers about 7% of the sun- 



lit surface area for a radius of 2 km, 45% for a radius of 
0.8 km, and rises to more than 80% for a radius of 0.60 km. 
Assuming that the nuclear radius of 46P/Wirtanen most 
likely is in the range 0.5-2.0 km we obtain from Eq. (2): 



[ .42 ± .01 .83 ± .01 1.66 ± .02 for Mod 1 
P rot [hrs] = I 3.0 ± 0.9 6.0 ± 1.8 12.0 ± 3.7 for Mod 2a 
[ 2.5 ±0.7 5.0 ±1.4 10.0 ±2.8 for Mod 2b 

where R e s is the effective radius of the nucleus. Thus, 
Model 1 gives very small rotational periods (fast rotation). 
Assuming that the mean density of the cometary material 
is equal to 1 g/cm 3 and the comet radius is equal to 1 km, 
one finds that the gravitational force at the nucleus sur- 
face is almost 20 times smaller than the centrifugal force 
that should cause rotational splitting of the nucleus. Alter- 
natively, Model 1 requires the comet radius to be greater 
than 4 km. 

Based on observations made at intermediate heliocen- 
tric distances and assuming a double-peaked light curve, 
(1997) found a possible rotational period 



Meech et al 



of about 7.6 hours. Lamy et al. (1998) derived a double 



peaked light curve with a rotational period of 6.0±0.3 hr. 
Models 2a and 2b are in agreement with these results. 
Assuming 6 hr for the rotational period our Models 2a 
and 2b give a radius in the range of 0.7-1.5 km and of 
0.9-1.7 km, respectively. 



From the range of variation of the light curve [Lamy et al 



(1998) found that the ratio of the axes must be greater 
than 1.20 for the assumed ellipsoidal body. Our values of 
1.14±0.05 (Model 2a) and 1.11±0.05 (Model 2b) for the 
R a /Rt ratio are only slightly below their lower limit. Thus, 
we concluded that Models 2a and 2b are to be preferred 
over Model 1 because of the smaller mean residual, and 
the reasonable prediction of the P rot /R a ratio. 

However, not all observational facts are explained by 
our models. The negative values of time shifts r (Model 
2b) and n, T2 (Model 2a) suggest that the gas produc- 
tion rate peaks before the perihelion passage. This is in 
contradiction with the positive light curve asymmetry 
given in Jorda & Rickman (1995) assuming that the jet 



force is produced only by sublimation directly from the 
nuclear surface. This question was also pointed out by 



Rickman & Jorda (1998). To explain this problem, they 



R eS [kr 



0.5 



1.0 



2.0 



discussed a model of 46P /Wirtanen in which the nucleus is 
subject to fragmentation along with sublimation. The icy 
chunks which break from the surface are expelled by the 
drag of the outflowing gas from the pores and propelled 
by a jet force due to sublimation. They argued that if the 
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outgassing from icy circumnuclear material is predomi- 
nant then it gives rise to the postpcrihclion brightness. 
Thus, it is possible that the gas production curve related 
to the nongravitational effects has no correspondence to 
the postpcrihclion shift of lightcurve maximum observed 
for 46P/Wirtanen. Furthermore, it is possible that the dis- 
continuity in global activity obtained from orbital calcu- 
lations would not be detected in the visible brightness of 
the comet. More studies are needed to judge whether the 
fragmentation hypothesis is true. 

For all of these reasons, we conclude that our forced 
precession model for the rotating comet nucleus (Model 
2b or 2a) reproduces - with very good precision - the 
nongravitational motion of 46P/Wirtanen during an al- 
most fifty-year interval and gives nucleus parameters (e.g. 
size and rotational period) consistent with the available 
photometric observations as well. Future evolution of the 
orbital elements and orientation of the comet spin-axis is 
presented in Table 2 for Model 2b since this model repre- 
sents the orbital motion with the smallest mean residual. 
One should note that the nongravitational forces have a 
significant influence on the cometary position even dur- 
ing a single apparition of the comet. Thus, these effects 
- even within a single apparition - must be included into 
any accurate orbit calculations. 
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